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Deep Chandra observations of TeV binaries - II: LS 5039 
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ABSTRACT 

We report on Chandra observations of the TeV emitting High Mass X-ray Binary 
LS 5039 , for a total exposure of ~70ks, using the ACIS-S camera in Continuos Clock- 
ing mode to search for a possible X-ray pulsar in this system. We did not find any 
periodic or quasi-periodic signal in the 0.3-0.4 and 0.75-0.9 orbital phases, and in a 
frequency range of 0.005 — 175 Hz. We derived an average pulsed fraction 3cr upper 
limit for the presence of a periodic signal of 15% (depending on the frequency and 
the energy band), the deepest limit ever reached for this object. If the X-ray emission 
of LS 5039 is due (at least in part) to a rotational powered pulsar, the latter is either 
spinning faster than ~ 5.6 ms, or having a beam pointing away from our line of sight, 
or contributing to ^ 15% of the total X-ray emission of the system in the orbital 
phases we observed. 
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1 INTRODUCTION 

The characterization of the compact objects hosted in the 
two TeV binaries LSI+61°303 and LS 5039 is now one 
long-standing question in high-energy astrophysics. The pe- 
culiar GeV (Abdo et al. 2009a) and TeV (Aharonian et al. 
2005, 2006) emission of these two High Mass X-ray Bina- 
ries (HMXBs) has been modeled as being powered either by 
the interaction of pulsar-accelerated particles in the pulsar 
wind zone (e.g., Sierpowska-Bartosik & Torres 2008) and/or 
by the interaction of shock-accelerated particles driven by 
the collision between the winds of the pulsar and the mas- 
sive companion star (e.g., Dubus 2006); or by a jet emitted 
by an accreting black-hole or neutron star (e.g., Khangulyan 
et al. 2008). No firm conclusion on this issue has been drawn 
thus far and the literature contains detailed models and dis- 
cussions on all of these alternatives. 

In Rea et al. (2010; Paper I) we have reported on a ~ 95 
ks Chandra observation of LS I +61°303 aimed at gathering 
the most stringent limit on (or detecting) the pulsed fraction 
of any X-ray signal from a putative energetic pulsar. In that 



paper we inferred an X-ray pulsed fraction for LS I +61°303 
of less than 10% (depending on the frequency and on the 
energy band), while the compact object was close to the 
apastron of its ~26 days orbit around the Be companion. We 
also observed two flares from LSI +61° 303 which showed a 
harder spectrum with respect to the source emission, and 
can be explained in an accretion scenario by variability in 
the accretion rate, while in the rotational-powered pulsar 
one, as the interaction between the pulsar wind and clumps 
in the Be wind. Similar flares have been observed by other 
experiments, among them, by RXTE (Smith et al. 2009, 
Torres et al. 2010, Li et al. 2011) and Swift-XRT (Esposito 
et al. 2006). 

The variation of the spectral index with flux has also been 
found for LS 5039 (Bosch-Ramon et al. 2005, Takahashi et 
al. 2009). However, the steadiness of the lightcurve is unique 
to this latter source, whereas there is significant orbit-to- 
orbit variability in LS I +61°303 (Torres et al. 2010, Li et al. 
2011). Kishishita et al. (2009) reported that LS5039 has a 



* Ramon y Cajal Research Fellow; rea@ieec.uab.es. 



2 Rea et al. 




0.2 



INT-2 



2xl0 4 3xl0 4 
Time (s) 



INT-3 




4xl0 4 




10* 1.5x10* 2x10* 2.5x10* 

Time (s) 



Figure 1. Chandra background-subtracted 0.3-10 keV lightcurves of the two segments of the Chandra observation of LS5039 binned 
at 100 s. 



X-ray lightcurve presenting long-term stability using multi- 
satellite observations between 1999 and 200o 

With these prior results we would not a priori expect 
to find in our Chandra observations anything different from 
this regarding the lightcurve evolution of LS 5039; our aims 
being rather the search for accretion and pulsation signa- 
tures. One possible explanation for the distinction in the 
lightcurve evolution of LS 5039 and LS I +61°303 is the fact 
that whereas the former system host an O star, LS I +61°303 
has a Be star primary, which is characterized with a dense 
equatorial disk. If, for instance, there are disc structure 
changes in time, the latter may induce variability in X-ray 
modulation, deviating its behavior for the clock-like emis- 
sion found for LS 5039 . 

The detection of fast X-ray pulsations would be an un- 
ambiguous tracer for a secure determination of the neutron 
star nature of the compact objects hosted in TeV binaries, if 
such is the case. Of course, in the accretion scenario Type I 
X-ray bursts might be another neutron star tracer, although 
never observed so far in a HMXB given the relatively high 
magnetic field of these young objects. In fact, for magnetic 
fields larger than ~ 10 10 G the thermo-nuclear instability 
giving rise to Type I X-ray bursts is suppressed, as the ac- 
creted H and He is burned stably in the accretion column 
of high B-field neutron stars (Lewin et al. 1995). Searches 
for pulsations from a young neutron star in massive binaries 
have in general many more chances of success in the X-ray 
than in the radio band. In fact, if the putative young neu- 
tron star is accreting, radio pulsed emission is expected to 
be halted by the accretion process, while if not accreting a) 
the X-ray pulsar beam is larger than the radio one, and b) 
the strong companion wind might well prevent the detection 
of radio pulsations because of free-free absorption and the 
large and highly variable Dispersion Measure (DM) induced 
by the wind (see e.g., Zdziarski et al. 2008). 



1 We note that Durant et al. (2011), using new Chandra obser- 
vations, report a flux ~2 times lower , and a smaller photon- 
index, than that observed by Suzaku (Kishishita et al. 2009) at 
the same orbital phase. This finding might represent a hint for 
a long-term variability, and they contradict the well-settled flux- 
gamma anti-correlation with the orbital phase (i.e. low flux always 
corresponded to large photon-indexes). We caveat that assuming 
the flux and spectrum of the Suzaku observation for the orbital 
phase covered by the Chandra one, we find that the pile-up in 
Chandra should be > 20%, and this might in principle account 
for the lower flux and harder photon index. 
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Figure 3. Search for the 3.5 ms spin period of PSR J0751+1807 
with PRESTO in two years of Fermi-hNT data. Top panel: here 
we did not correct the ToAs for the orbital period, and the pulsar 
spin period cannot be detected. Bottom panel: here we correct the 
ToAs for the orbital period, detecting a clear signal at the pulsar 
rotational period. 



Deep searches for pulsations have been performed in the 
radio band at several frequencies, with the hope of detecting 
a fast spinning radio pulsar as in the case of the other TeV bi- 
nary, PSRB1259-63 (Johnston et al. 1999, 2005). However, 
no radio pulsation has been detected so far from any of these 
two sources. In that respect it is crucial to note that at pe- 
riastron PSR B 1259-63 does not show radio pulsations. Its 
periastron (given the large orbit, 3.4 year period) has about 
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Figure 2. Left panel: LS 5039 's geometry considering the orbital solution of Aragona et al. (2009), and the phases for Inferior conjunction 
(INFC), Superior conjunction (SUPC), periastron (P), and apastron (A) are marked accordingly, but the inclination is not taken into 
account in the plot. The orbit (black solid line) and the massive star (MS; in orange) are roughly to scale. The blue thick line indicates 
the orbital phases spanned by our two Chandra observations. Right panel: Chandra data folded with the orbital ephemeris of LS 5039 . 



the same dimension of the major axis of LS I +61°303's or- 
bit (26 days period) and it is way larger than LS 5039 's obit 
(4 days period), hence it is reasonable to expect that no ra- 
dio pulsed emission is observed from these much compact 
binaries. 

Before the Chandra observations we report here (and in 
Paper I for LSI+61°303), archival observations which could 
give reliable upper limits on pulsations for a possible fast 
spinning pulsars (P ^ 100 ms) hosted in LSI+61°303 and 
LS 5039 were not very constraining, coming mainly from 
RXTE and XMM-Newton observations. In particular the 
high background of these instruments (especially .RXTE) 
limited the pulsed fraction sensitivity of these observations. 
For LSI +61° 303 the deepest pulsed fraction limit was de- 
rived from a 41 ks XMM-Newton observation (Sidoli et al. 
2006), which was < 28%, in the 12-200 ms period range. 
Note that the RXTE monitoring observations performed in 
1996 gave an upper limit of only 32% in the 1-200 ms range, 
in fact, even though RXTE's timing resolution and collecting 
area was much larger than XMM-Newton, the much higher 
background contamination is limiting the sensitivity to weak 
signals, resulting in a larger pulsed fraction limhff]. Similarly, 
for LS 5039 the deepest limits for the presence of a fast 
pulsar were derived from a 50 ks RXTE observation per- 
formed in 2003, giving an upper limit of the pulsed fraction 
of < 30%. 

We present here the Paper II of our series of Chandra 
observations aimed at searching for pulsed X-ray emission 
from TeV binaries. In Paper I we presented such analysis 
for a ~95ks Chandra Continuous Clocking mode observa- 

2 For that observations, Harrison et al. (2000) claim a limiting 
pulsed fraction of ~6%. However they considered the total count 
rate without correcting for the cosmic and instrumental back- 
ground, which if corrected increases substantially the upper limit 
on the detectable pulsed fraction. 



tion of LS I +61° 303, while in this paper we present a similar 
analysis on ~70 ks of Chandra data on LS 5039 . In all obser- 
vations we used the same instrumental set-up, to favor the 
search for fast spinning pulsars. In §2 we report on the ob- 
servation and analysis, we present our results in §3, followed 
by a discussion (§4). 



2 OBSERVATION AND DATA ANALYSIS 

The Advanced CCD Imaging Spectrometer (ACIS) camera 
on board of Chandra observed LS 5039 twice, on 2009 July 
31st (start time 05:41:49 UT; Obs-ID 10053) and on 2009 
August 6th (start time 06:23:48; Obs-ID 10932) for an ex- 
posure time of 45.138 and 25.079 ks, respectively. The obser- 
vations were performed in Continuos Clocking (CC) mode 
(FAINT), providing a time resolution of 2.85 ms and imag- 
ing along a single direction. In both pointings the source 
was positioned in the back-illuminated ACIS-S3 CCD at the 
nominal target position. Standard processing of the data was 
performed by the Chandra X-ray Center to level 1 and level 
2 (processing software DS ver. 8.0). The data were repro- 
cessed using the CIAO software (ver. 4.2) and the Chandra 
calibration files (CALDB ver. 4.2). 

We corrected the times for the variable delay due to 
the spacecraft dithering and telescope flexure, starting from 
level 1 data and assuming that all photons were originally 
detected at the target position. Furthermore, data were fil- 
tered to exclude hot pixels, bad columns, and possible after- 
glow events (residual charge from the interaction of a cosmic 
ray in the CCD). Photon arrival times are in TDB and were 
referred to the barycenter of the Solar System using the 
JPL-DE405 ephemeris. 

In order to carry out a timing analysis, we extracted 
the events in the 0.3-10 keV energy range from a region of 
5x5 arcseconds around the source position (RA 18:26:15.00, 
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Figure 4. Left panels: Top row: power spectra of LS 5039 in the 0.3-10 keV energy band binned at 0.00285 s with superimposed the 3<r 
upper limits on the detection of a periodic signal on top of the Poisson noise (we used here logarithmic scales to better show the low 
frequency range). Bottom row: pulsed fraction 3cr upper limits on the semi-amplitude of a sinusoidal signal. Right panels: pulsed fraction 
3<r upper limits versus frequency in three different energy ranges (from top to bottom: 0.3— 2keV, 2— 4keV, 4— lOkeV), calculated on the 
power spectrum binned at .00285 s (left panel) in the 0.005—175 Hz frequency range (see H3.2l for further details). 



Dec —14:50:53.6). The source spectrum was extracted from 
a rectangular region of 5x25 arcseconds around the source 
position, and the background was taken independently from 
a source-free region in the same chip. We extracted the 
appropriate response matrix files (RMFs) and ancillary re- 
sponse files (ARFs) for each observation. The source ACIS- 
S average count rate in the 0.3-10 keV energy band was 
0.509 ±0.003 counts s _1 and 0.430 ± 0.004 counts s _1 for the 
first and second observation, respectively (see Fig.[T}. 



3 RESULTS 

Using the orbital ephemeris of Aragona et al. (2009), namely 
an orbital period of 3.90608±0.00010 days and the zero- 
phase assumed at T =52825.485±0.053 MJD, we calculated 
that the first Chandra observation was performed while the 
compact object was passing from phase 0.3 to 0.4, while 
the second observation was spanning phases 0.75-0.9 (see 
Fig.®. 

3.1 Orbital demodulation 

Given the relatively long Chandra exposure time with re- 
spect to the ~ 4 days orbital period of our target, we had 
to correct the photon arrival times for the orbital motion of 



the system, before searching for any periodic signal (which 
might otherwise be washed out; see e.g. Vaughan et al. 1994 
and reference therein). To correct the arrival times of each 
X-ray photon for the orbital motion of LS 5039 we devel- 
oped an algorithm in the framework of the program TEMP02 
(Hobbs et al. 2006). Unfortunately, a dedicated plug-in for 
the orbit demodulation is not present in the official release, 
hence we built our own program for this purpose. In par- 
ticular, we used the TEMP02 output format general2, that 
calculates the pulsar pulse number and phase corresponding 
to each arrival time, giving as input the orbital ephemeris 
of a binary system. The algorithm we developed is based on 
the creation of a fake TEMP02 parameter file that includes the 
orbital parameters of LS 5039 (Aragona et al. 2009), and a 
pulsar with a constant spin frequency, which we assumed 
being 1 Hz. In this way the 1 Hz pulsar works like a clock on 
the reference frame comoving with the true compact object, 
and the values returned by the output format general2 cor- 
respond to the photon arrival times in the pulsar reference 
frame, hence demodulated by the orbital period. 

We have tested this algorithm both on simulated and 
real data. We have first simulated with the TEMP02 plug-in 
fake a 330 ms pulsar in a binary system with 4 days or- 
bital period. In particular, we simulated the Time Of Ar- 
rivals (TOAs) of the pulsar with a Gaussian random noise 
with a root mean square (r.m.s.) set to 0.1ms. We then 
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Figure 5. Left panel: spectra of the 5 time intervals fitted with an absorbed (Nh = (7.1 zt 0.4) X 10 21 cm — 2 ) power-law function (see 
FigfjJ and text for details). Right panel: la sensitivity for the presence of spectral lines (this corresponds to the limit on the flux of 
detectable spectral lines as broad as the energy-dependent instrument resolution on top of the continuum flux as a function of energy) 
for each for the five spectra. 



checked that the residual plot of the simulated TOAs using 
the full orbital ephemeris we have input was identical to the 
residual plot of the TOAs derived demodulating the orbital 
motion using the above algorithm, when inputting only the 
pulsar parameters. As a further check, we have applied our 
demodulation algorithm on two years of Fermi Large Area 
7-ray Telescope (LAT) data of the binary millisecond pul- 
sar (PSRJ0751+1807; Nice et al. 2005; Abdo et al. 2009b), 
which is similarly composed by a 3.5 ms pulsar in a 0.26 days 
orbital period with an helium-rich white dwarf. In Fig. [3] we 
show the PRESTO prepf old outcome of the Fermi-LAT data 
before (top) and after (bottom) applying the demodulation 
algorithm, showing that after correcting for the orbital pa- 
rameters we can detected the spin period of the pulsar. 

For the orbital parameters we use: an orbital 
period of 3.90608±0.00010 days, a projected semi- 
major axis of 3.329±0.153 s, a periastron longitude of 
236.0±5.8 degrees, and an eccentricity of 0.337±0.036, 
referred to T =52825.485±0.053 MJD0. 

3.2 Pulsation search and light-curve analysis 

Using the event files with the photon arrival times referred 
to the Solar System barycenter and corrected for the or- 
bital motion (see above), we searched for periodic and quasi 
periodic signals with unprecedented sensitivity, thanks to 
the extremely low ACIS-S background. Given the length 
of our observations (70 ks in total), the timing resolution 
of the CC-mode (0.00285 s), and the number of counts of 
our observation, we could search for periodic signals in the 
~ 0.005 — 175 Hz frequency range. As we explained in detail 
in Paper I, we studied the source power spectra performing 

3 For technical reasons, these parameters in Aragona et al. (2009) 
are referred to an epoch expressed in Heliocentric Julian Day, 
that differs from the most commonly used Julian Days by a vari- 
able amount of time with a maximum difference of 500 s. Given 
~5000 s error in the epoch determination we neglected the differ- 
ence between the Heliocentric JD and the JD we assume in this 
work. 



Fast Fourier Transforms (FFTs), considering the total ex- 
posure time (summing both observations), and on each of 
the five time intervals shown in Fig.[T] Furthermore, we per- 
formed the same search dividing the entire data set in three 
energy bands (0.3-2, 2-4, and 4-10 keV). 

For computing reasons, we performed an average over 
13 FFTs with a bin time of 0.00285 s (Fig.g] top-left panel), 
resulting in 2,097,152 bins (~ 6ks) for each of the 13 av- 
eraged power spectra (the averaged power spectrum had 
a xt distribution with 26 degree of freedom (d.o.f.)). For 
the power spectra produced for the three different energy 
bands, given the lower number of counts, we could use a 
single FFT (their power spectrum had then a xS distribu- 
tion with 2 d.o.f.) . We took into account the number of 
bins searched, and the different d.o.f. of the noise power 
distribution (Vaughan et al. 1994; Israel & Stella 1996) in 
calculating the 3<r detection upper limits reported in Fig.g] 
. We accounted for the red-noise (despite negligible in this 
observation) in the calculation of the detection significance 
by using a smoothing window technique (see Israel & Stella 
1996 for further details). We did not find any periodic or 
quasi-periodic signal. 

As reported in Paper I, we computed the 3<r upper lim- 
its on the sinusoid amplitude pulsed fraction (PF), com- 
puted according to Vaughan et al. (1994) and Israel & Stella 
(1996), which ranges in the 0.3-10 keV energy band between 
PF <13-21% (see Fig.g] bottom- left panel). In Fig.g] (right 
column) we calculate the same limits as a function of the 
energy band, which given the lower number of counts causes 
the energy-dependent PF limits to be slightly larger than 
those derived using the whole datasets. 

In addition, we looked for the presence of SGR-like 
bursts by binning the counts in intervals of 0.1 s and search- 
ing for excesses above a count threshold corresponding to 
a chance occurrence of 0.1% (taking into account the total 
number of bins), but we did not find any significant short 
burst. 

In Fig.[T] and [2] we report on the light-curve of both 
observations, which show a variability correlated with the 
orbital phase of the system, as previously observed in a long 
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Figure 6. Le/4 panel: Pulsed fractions versus total X-ray luminosity. Right panel: Pulsed fractions versus rotational power. See also Tab. 1. 
Circles stand for sources with associated TeV conterparts. Down arrows and horizontal lines refer to the upper limits for LSI+61°303 
and LS5039 . See §4.2 for details. 



Suzaku observation (Kishishita et al. 2009; Takahashi et al. 
2009). 

3.3 Spectral analysis: continuum modeling and 
search for spectral lines 

All spectra reported here were binned such to have at least 
25 counts per bin and not to over-sample the instrument 
spectral resolution by more than a factor of 3. We fit the 
spectra relative to the first and to the second observation 
with an absorbed power-law, which give a xt = 1-18 an d 
xi = 1-32 for 421 and 261 d.o.f., respectively. Given the 
known spectral variability of LS 5039 as function of orbital 
phase, we then searched for a spectral evolution dividing the 
observation in 5 time intervals (see Fig. [T] and [S]). Fitting to- 
gether the spectra of the 5 time-intervals with an absorbed 
power-law function we found a xt = 1-15 for 725 d.o.f, with 
the absorption value tied to be the same for all 5 spectra and 
the slope and normalization of the power-laws free to vary. 
We find an absorption value of N H = (7.1 ±0.4) x 10 21 cm~ 2 
(with abundances from Anders & Grevesse (1989), and 
cross-setion for photoelectric absorption from Baluchinska- 
Curch & McCammon (1998)), and a power-law index which 
varies only very slightly, from 1.49 to 1.55 depending on the 
time-interval (see Fig. [5] all errors in the spectral values 
are reported at la confidence level). The source 0.5-10 keV 
absorbed flux varies from the first to the fifth time interval 
(from left to right in Fig.[TJ: 8.8 ± 0.2, 7.8 ± 0.2, 6.7 ± 0.2, 
6.0±0.2 and 7.1±0.2x 10~ 12 ergs _1 cm~ 2 (corresponding to 
unabsorbed fluxes of 1.2±0.2, 1.0±0.2, 0.9±0.2, 0.8±0.2 and 
1.0±0.2 x 10 _11 ergs _1 cm -2 ). Furthermore, we searched for 
possible evidence for Nh variability among the five spectra 
by means of contour plots with respect to the photon index 
(r), but while the spectral indexes where found to show a 
subtle variability (at 2.5cr in the best case, between interval 
3 and 5; see also Fig.JSJ, we do not have enough counts to 
constrain any variability in the absorption value. 

Comparing our flux and spectral parameters with those 
derived by Suzaku at similar orbital phases (Kishishita et al. 
2009; Takahashi et al. 2009), we find a perfect match within 
a la error (see Tab. 2 in Takahashi et al. 2009). 



A deep inspection of all spectra did not show evidence 
for significant absorption or emission features, for which we 
report in Fig.[5] (right panel) the la limit on the flux of 
detectable spectral lines as broad as the energy-dependent 
instrument resolution on top of the continuum flux, as a 
function of energy (we refer to the XSPEC modeling package 
for detailsQ). 



3.4 Spatial analysis: search for extended emission 

An ~1' extended X-ray emission around LS 5039 has been 
recently reported by Durant et al. (2011). We produced 
a one-dimensional radial profile of LS 5039 for both the 
Chandra observations, using the CC mode one-dimensional 
imaging capability, to search for any evidence for such dif- 
fuse extended X-ray emission. We first generated an image 
of the one-dimensional strip in the 0.3-10 keV band, and 
then subtracted the background count rate to remove in- 
strumental and cosmic X-ray background. We extracted a 
one-dimensional surface brightness distribution extracting 
photons from 125 annuli centered on the source position 
and of 2 pixel wide (~0.99"). We then did a similar extrac- 
tion for the simulated Chart/MARX Point Spread Function 
(PSF). We did not find evidence for a significant extended 
emission or structures within ~1' around our target. 

Assuming an absorbed power law spectrum with Nh = 
6.4 x 10 21 cm~ 2 and T = 1.9 for the extended emission be- 
tween 20-60" (Durant et al. 2011), we inferred a 3<r limiting 
unabsorbed flux for the presence of X-ray diffuse emission 
of ~ 3 x 10 -13 in the 0.5-10 keV energy range. The limit we 
find is higher than the flux of the X-ray extended emission 
reported for LS5039 (1.4 x lO'^ergs^cm" 2 ; see Tab 2 in 
Durant et al. 2011), hence it is compatible with the presence 
of such emission, and we cannot constrain any further this 
emission with the observations reported here. 



4 http://heasarc.nasa.gov/xanadu/xspec/ 
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4 DISCUSSION 

We reported here on a ~70ks Chandra look at the TeV bi- 
nary LS 5039 aimed at searching for a dim periodicity in its 
X-ray emission due to the presence of a pulsar in the sys- 
tem. As in the case of LSI+61°303 (Paper I), we did not 
find any periodic or non-periodic signal in these observa- 
tions, which span a significant fraction of the orbit (namely 
phase 0.3-0.4 and 0.75-0.9). We confirm the orbit-by-orbit 
stability of the light-curve shape and the modest spectral 
variations previously observed in this system (Bosch- Ramon 
et al. 2005; Kishishita et al. 2009; Takahashi et al. 2009). We 
derived upper limits on the X-ray pulsed fraction of LS 5039 
(PF 15%; see §3.2\ and Fig.[4|, as a function of frequency 
and energy band. 

4.1 On the impact of different ephemeris for 
LS 5039 

A recent analysis of the mass and orbital parameters of 
LS5039 was presented by Sarty et al. (2010), using obser- 
vations obtained with the MOST satellite. These authors 
report two main findings: 1) the orbit of LS 5039 has a sig- 
nificantly lower eccentricity than reported before, and 2) the 
total mass of the binary system is toward the higher end 
of previous estimates. A direct comparison of their results 
with respect to those of Aragona et al. (2009) and Casares 
et al. (2005) shows that despite the lower mean value of the 
eccentricity found by Sarty et al. (2010), all of them are con- 
sistent within la error (see Table 1 in Sarty et al. 2010). All 
of the other orbital parameters are in agreement with those 
of Aragona et al. (2009) within errors, except for the system 
mass function, the value of which is instead compatible with 
that obtained by Casares et al. (20050- For completeness we 
performed our timing analysis (see § 13. H and SI3.2I) using the 
three ephemeris, finding similar results. 

4.2 Comparison of our limits with X-ray pulsed 
fractions of known pulsars 

We use the deep PF upper limits we derived here for 
LS5039 (see andFig.[i, and in Paper I for LS I +61°303, 
to draw a comparison between the X-ray emission of the pu- 
tative pulsar in these TeV binary systems, with that of ac- 
creting pulsars in HMXBs, and isolated rotational-powered 
pulsars. In the HMXB case, after a thorough search in the 
literature, the large variability of such pulsed fractions (from 
a few % to 90%) as a function of luminosity and energy even 
in a single binary system, did not allow us to draw any con- 
fident conclusion from the comparison with our limits. The 
large variability of the PF in these systems is probably due 
to the variable accretion rate, and accretion column geome- 
try. 

On the other hand, in Tab. 14. II and Fig. [5] we report on 
the comparison between the X-ray PF values of all isolated 
rotational-powered pulsars where this value has been mea- 
sured, and the limits we have for LS I +61° 303 and LS 5039 . 

5 f(m) (Mg) =0.0053 ± 0.0009, 0.0026 ± 0.0004, 0.0049 ± 0.0006 
for Casares et al. (2005), Aragona et al. (2009), and Sarty et al. 
(2010), respectively. 



All pulsed fractions here are calculated as: 

where N max and N m in are the maximum and minimum 
number of counts in the X-ray folded light-curve, respec- 
tively. We have tried to use the same X-ray band for all 
sources, although not for all of them the PF value was re- 
ported in the 0.3-10 keV energy range we wanted. In these 
few cases where only the 0.3-2 keV or 2-20 keV range were 
available (as e.g. for J1930+1852 which was observed by 
itXTE), we have checked that the PF values had only a 
marginal dependence with energy, and assumed that the 
same value did hold also for the 0.3-10keV energy range. 
Furthermore, for those sources for which in the original pa- 
per the PF was given using a different definition, we have 
extracted the values from the pulse profile directly, such 
to compare PF values derived as coherently as possible. In 
Fig. [6] we plot the values reported in Tab. 14. II using red cir- 
cles for pulsars being associated with a TeV source (see also 
Mattana et al. 2009), and black triangles for pulsars with 
no detected TeV counterpart. Note that the TeV emission 
is mainly related to the presence of a pulsar wind nebula 
sustained by the energetic pulsar. 

The luminosities we report in Fig. |S] for LSI+61°303 
and LS5039 are 7.6 and 7.5 x 10 33 ergs _1 (they assume 
a distance of 2.3 kpc and 2.5 kpc, respectively), and were 
derived from the average source flux reported in this work 
and Rea et al. (2010, Paper I), hence they refer to the orbital 
phase ranges where such observations were performed. 

At variance with what was stated in Paper I, although 
it is true that isolated rotational powered pulsars have on 
average PF larger than the limits we found for TeV binaries, 
the current detailed analysis on their X-ray pulsed fractions 
using a coherent definition of ^" tax ~^" li " for all pulsars, 
does not seem to show a trend of higher PF values for sources 
with TeV counterparts. 
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Name 


log E rot 


log L x tot 


B-field 


Period 


PF 


Reference 




erg/s 


erg/s 


G 


ms 







B0531+21 (Crab) 


38.65 


37.04 


3.78e+12 


33.40 


99.3% 


Becker ct al. 2002; Tennant et al. 2001; Li et al. 2008 


B0833-45 (Vela) 


36.84 


32.83 


3.38e+12 


89.29 


6.6% 


Becker ct al. 2002; Halpern ct al. 2001; Li ct al. 2008 


B0633+17 (Geminga) 


34.51 


31.10 


1.63e+12 


O O T / 1 M 

237.09 


40% 


Becker ct al. 2002; Caraveo et al. 2004 


B1509-58 (MSH 15-52) 


37.25 


34.29 


1.54e+13 


150.23 


67% 


Becker et al. 2002; Kuster et al. 2002; Li et al. 2008 


B1929+10 


33.59 


30.00 


5.18e+ll 


226.51 


39.0% 


Becker ct al. 2002; Slowikowska ct al. 2005 


PSR B0950+08 


32.75 


29.92 


2.44e+ll 


253.1 


33.0% 


Zavlin ct al. 2004 


B1821-24 


36.35 


33.24 


2.25c+09 


3.05 


100% 


Becker et al. 2002; Rutlcdgc ct al. 2004 


B0656+14 


34.58 


32.98 


4.66e+12 


O O A O T 

384.87 


29% 


Becker et al. 1997; Chang et al. 2001 


B0540-69 


38.17 


36.39 


4.98e+12 


50.37 


50% 


Becker ct al. 1997; Kaarct ct al. 2001 


TOT O A OO 

J2124-33 


33.55 


30.35 


00 1 r\o 

3.22c+08 


4.93 


177o 


Becker ct al. 1997; Bogdanov et al. 2008 


B1055-52 


34.48 


33.42 


1.09e+12 


197.10 


25% 


Becker et al. 2002; Dc Luca ct al. 2005 


T AO I 1 A O'lO 

JU21o+42oz 


35.37 


32.75 


a on — 1 no 


2.32 


63% 


Becker ct al. 2002; Kuiper ct al. 2002; Mineo et al. 2000 


PSR J 1617-5055 


37.20 


33.26 


3.10e+12 


69.33 


47% 


Becker ct al. 2002; Kargaltsev et al. 2009 


PbK J0030+0451 


33.53 


30.56 


tr „ 1 no 

2.25C+08 


4.86 


337o 


Becker et al. 2002 


B1937+21 


36.04 


31.66 


4.09e+08 


1.55 


92% 


Becker et al. 2002; Nicastro et al. 2002 


J 0205+6449 


37.42 


32.31 


3.61e+12 


65.67 


72% 


Becker ct al. 2002: Murray ct al. 2002 


J2229+6114 


37.34 


33.25 


2.03e+12 


51.62 


197o 


Becker ct al. 2002; Halpern ct al. 2001 


J1930+1852 


37.08 


34.53 


1.03e+13 


136.85 


2.5% 


Becker ct al. 2002; Leahy ct al. 2008a; Lu ct al. 2007 


J1811-1926 


36.80 


30.29 


1.71c+12 


64.67 


26% 


Becker ct al. 2002; Torii ct al. 1999 


J0537-6909 


38.68 


36.30 


9.25e+ll 


16.11 


48% 


Becker ct al. 2002; Marshall ct al. 1998, Wang et al. 2001 


J1420-6048 


37.00 


29.35 


2.41e+12 


68.18 


35% 


Becker et al. 2002; Roberta et al. 2001 


PSR J1838-0655 


36.70 


30.92 


1.89e+12 


70.5 


4.8% 


Gotthelf et al. 2008; Lin et al. 2009 


PSR ,11846-0258^ 


36.9 


34.89 


4.86e+13 


325.7 


14% 


Kuiper ct al. 2009; Leahy ct al. 2008b 


PSR J1119-6127 - 


36.40 


33.00 


4.10e+13 


407.7 


51% 


Gonzalez et al. 2005 


B0355+54 


34.65 


30.92 


8.39e+ll 


156.38 


43% 


McGowan et al. 2007 


B0628-28 


32.17 


30.00 


3.01e+12 


1244.42 


50% 


Tepedelenhoglu et al. 2005; Li et al. 2008 


J1124-5916 


37.07 


34.53 


1.02e+13 


135.31 


28% 


Hughes ct al. 2004; Li ct al. 2008 


B 1706-44 


36.53 


32.65 


3.12e+12 


102.46 


30% 


McGowan ct al. 2004; Li ct al. 2008 


J1747-2958 


36.39 


34.38 


2.49e+12 


98.81 


29% 


Hickox ct al. 2009; Li ct al. 2008 


B1951+32 


36.56 


33.54 


4.86e+ll 


39.53 


1.8% 


Chang et al. 1997; Li et al. 2008 


J1357-6429 


36.49 


32.30 


7.83e+12 


166.1 


73% 


Zavlin ct al. 2007 



Table 1. Parameters for all isolated rotational-powered pulsars with detected X-ray pulsations. The _E ro t is the rotational power as 
derived from the rotational period and its derivative; L x tot is the total X-ray luminosity, and the PF is defined here as j^ max ~j^""" , 
where N ma:E and N m i n are the maximum and minimum number of counts in the folded light-curve. We refer to the quoted papers 
for more details on L x tot and PF, which were always in an energy range close to 0.3-10keV. The values of the the rotational power 
(Erot), B and the period were all derived from the ATNF pulsar catalog (Manchester et al. 2005). See also Fig. [6] 



a Note that PSR J1846-0258 has shown magnetar-like activity (Gavriil et al. 2008; Kumar & Safi-Harb 2008), hence it should also 
have magnetic power contributing to its X-ray emission, rather than being "only" rotational powered. 
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